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I. INTRODUCTION 

This is a sia - month iirogress report covering research on the moon and 
its environment by the principal Investigator. The specific objectives, ^ 

approaches, and statements of work are presented in the body of the proposal 
submitted in February 1975* Section II briefly summarizes results in the 
primary research areas, and Section III lists publications and verbal papers 
submitted during the six-month period. Copies of publications and abstracts 
are included in Section IV. 

II. SUMMARY OF RESULTS 

During the time period 1 September 1975 - 28 February 1976 our lunar 
analysis has expanded to include data from a total of nine lunar magnetome- 
ters: three Apollo lunar surface magnetometers, two Apollo subsatellite 
magnetometers, two magnetometers aboard Explorer 35^ and the Russian Luna-22 
orbiter and Lunokhod-2 surface rover. New results have been obtained concern- 
ing electromagnetic, compositional, and structural properties of the lunar 
interior. In particular, the following tasks were con^pleted during the grant 
period. 

(l) Lunar Magnetic Permeability and Iron Abundance : Recent research 

studying magnetic permeability and iron abundance of the moon has involved 
the use of two different analytical techniques and a total of five different 
instruments. In the first technique simultaneous data are used from one 
Apollo surface magnetometer and one or both magnetometers aboard the lunar 
orbiting Explorer 35* A total of 10 lunations of quiet geomagnetic tail 
data have been used, extending over parts of 3 different calendar years. 

In the second (more recent) technique we use 5 lunations of simultaneous 
data from two Apollo surface magnetometers. This technique uses a differen- 
tial method of subtracting out the external magnetizing field and has tne 
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advantages of (a) using all three vector components in the analysis rather 
than just the radial component, (b) subtracting out all constant fields 
measured out at either site, mailing its anailysis independent of offset 
errors and remanent field determination at either site, and (c) making use 
of hl^ resolution LSM data only. The table below summarizes magnetic per- 
meability results to date. 


Magnetometers 

Wo. Lunations 

Bate Components 

Used 

Permeability 

AP 12 LSM/EX 35 

8 


1. 012±0. 006 

AP 15 I£M/EX 35 

2 


1.019±0.004 

AP 15 LSM/aP 16 LSM 

5 

1972 all 

1.008+0.005 


We are able to calculate iron ab\mdance in the moon from the bulk 
lunar permeability by assuming a suitable mineral compositional model of 
the moon (olivine or orthopyroxene of density and percentage iron silicate 
composition consistent with known global density, moment of inertia & other 
geochemical constraints); and a Curie isotherm depth consistent with lunar 
temperature pi-ofiles determined by seismic, magnetometer, jmd geochemical 
results. The free iron content is found to be 2.5 ± 2.0 percent by weight 
and the total iron in the moon is 9*0 ± 4.7 wt 

(2) Limits on a Highly-Conducting Lunar Gore ; The moon could have a 
small highly conducting Ivtnar core which, if having a sufficiently long 
Cowling time of the order of a few days or longer, would continually exclude 
the external geomagnetic tail field, therefore acting as a "diamagnetic" 
region of effective permeability =0. In this case the field enhancement 
of the outer region dominated by ferromagnetic and paramagnetic iron would 
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compete with the field decrease of the "diamagnetic” corej magnetometers on 
the surface or in lunar orbit would then measure the resultant net field. 

Our previously reported results assumed no such diamagnetic core; indeed, 
nothing in our permeability or conductivity studies indicates the existence 
of such a core. 

A preliminary calculation of the upper limits on the size of a highly 
conducting (assumed 0 = “) core has been completed for an assumed Curie iso- 
therm depth of 250 km (consistent with electrical conductivity results). We 
have calculated two cases for lunar compositions of olivine and orthopyroxene, 
in each case assuming no iron silicate, so that free iron content is maximized 
and therefore size of the hypothetical core is maximized. The upper limit 
on an orthopyroxene moon is 535 km and for an olivine moon 385 hm. It is 
noted that for both lunar models the minimum core radius is zero, that is, 
qii our measurements to date are consistent with the nonexistence of a highly 
conducting core. 

(3) Lunar Electrical Conductivity ; A new analytical technique has 
been recently developed and applied to lunar electrical conductivity, in 
which simultaneous data measured in high-latitude regions of the geomagnetic 
tail are used from a network of three instruments: the Apollo I 5 lunar sur- 

face magnetometer (IfiM), the Apollo I 6 LSM, and the Apollo I 6 subsate3JLite 
magnetometer, which provides coverage around the entire global circumference. 

A numerical analysis program uses measured total fields from the three mag- 
netometers to separate the time-varying field from the induced lunar poloidal 
eddy current field. Thereafter a conductivity profile for the lunar interior 
can be found by an iterative technique. 
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Also^ a two “instrument teclmique has been applied whicli uses Apollo I'd 
I^M data and simultaneous data from the Goddai’d magnetometer aboard Explorer 
35. All measurements are made when the moon is located in high-latitude 
regions of the geomagnetic tail where plasma effects in the lunai’ environ- 
ment ari minimal. Individual magnetic events are superimposed to obtain 
a single large transient for analysis. Also, during examination of five 
years of data one exceptionally large, well-behaved transient was recorded 
when the moon was in the geomagnetic tail. This single event promises sub- 
stantial improvement in resolution and sounding depth for conductivity 
analysis. Most recent electrical conductivity results are included in 
Section IV. 

(4) Lunar Internal Structure Inferx*ed from Conductivity and Permeability 
Results ; Structural implications of conductivity results to date (Dyal et 
al. , 1976) are as follows: The first 250 Ion depth has very low electrical 
conductivity. No conductivity transition is seen at the 60 km seismic dis- 
continuity reported by Nakamura, Latham and coworkers, but our resolution 
to date is limited at such shallow depths. Between 250-300 km depth the con- 
ductivity increases dramatically, reaching a "knee" at about 300 km depth. 

This change corresponds closely to the seismic discontinuity at that depth. 

From 300-900 km depth the conductivity continues increasing, but at a 
decreased rate. At about 900 km a slight increase in conductivity is indicated 
by the best fit profile at the seismic velocity change in the moon quake region 
of the lunar interior, but errors on the conductivity profile are large enough 
at these depths to prevent a definite conclusion at this time. A highly con- 
ducting core of maximum radius 535 hm is found to be possible for the extreme 
case of magnesium-silicate dominated orthopyroxene with Curie isotherm depth 
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of 250 kin indicated from the temperature profiles inferred from conductivity 
calculations. However^ there is no positive indication at this time that 
any core of conductivity > 10 mho/m need exist in the moon. 

Ill, LIST OF PUBLICATIONS /abstracts 

(1) I)yal P. > Parkin C.W. , and Daily W.D. , "Lunar Electrical Con- 
ductivity and Magnetic Permeability.” Proc. Lunar Sci. Conf . j 
6thj vol. 3, pp. 2909-2926, 1975* Pergamon Press. 

(2) Dyal P. , Parkin C.W. , and Daily W.D. , "Structure of the Lunar 
Interior from Magnetic Field Measurements," in Lunar Science VII , 
part I, 1976. 

(3) Daily W.D., Dyal P. , and Parkin C.W. , "Lunar Electrical Conducti- 
vity Profile from Magnetometer Network," (abs.) to be published 
in EOS , vol. 1976. 
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Ivor. Lunar S(i. Conf. m :9CW0926. 

Printed in the Uniied Sintci of Amcrici 


Lunar electrical conductivity and magnetic permeability 

Palmer Dyal 

NASA Ames Research Center, Motlctt Field, California ^4035 

Curtis W. Parkin 

Department of Physics, University of Santa Clara, Santa Clara, California 950S3 

William D. Daily 

Department of Physics and Astronomy. Brigham Young University, Provo, Utah 84602 

Abstract— Improved analytical techniques are applied to a larger Apollo magnetometer data set to 
yield values of electrical conductivity, temperature, magnetic permeability, and iron abundance. 
Average bulk electrical conductivity of the moon is calculated to be 7 x 10'' mho/m. Allowable 
solutions for electrical conductivity indicate a rapid increase with depth to — t0'’mhu/m within 
230 km. The upper limit on the calculated size of a hypothetical highly conducting core |a>7x 
10*’ mho/m) is O.SIRIR^^. The temperature profile, obtained from the electrical conductivity profile 
using the laboratory data of Duba et at. (1974) for olivine, indicates high lunar temperatures at 
relatively shallow depths. These results imply that the Curie isotherm is at a depth of less than 200 km. 
Magnetic permeability of the moon relative to its environment is calculated to be 1.008 ±0.005, 
Adjustment of this r»$uli to account for a diamagnetic lunar ionosphere yields a lunar permeability, 
relative to free space, of 1.012 traa. Lunar iron abundances corresponding to this permeability value 
are 2.5 2’:’ wt.?c free iron, and 5.0-13.5 wt.% total iron for a moon composed of a combination of free 
iron, olivine, and orthopyroxene. 


Introduction 

Data from the network of magnetometers placed on the moon by Apollo 
astronauts have allowed investigation of internal lunar properties. The purpose of 
this paper is to report our latest results of lunar electrical conductivity, tempera- 
ture, magnetic permeability, and iron abundance. 

Previous electrical conductivity analyses using a time-dependent transient 
response technique have been applied to nightside lunar data with the moon in the 
solar wind (e.g., Dyal and Parkin. I97ta; Dyal et al„ 1972), dayside data in the 
solar wind tDyal et ai, 1973). and to geomagnetic tail data (Dyal et aL 1974), 
Results to be presented here are calculated using data from deep in the 
geomagnetic tail lobes. A new transient-superposition technique is used which 
improves the signal-to-noise ratio in the analysis. In addition, amplitude and phase 
information are used for all three vector components of magnetic field data. 

In our earlier work on magnetic permeability and iron abundance of the moon 
(Dyal and Parkin, 1971b; Parkin et at., 1973, 1974), the analytical technique 
involved use of simultaneous data from the lunar orbiting Explorer 35 mag- 
netometer and the Apollo 12 or IS surface magnetometer. With the present 
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technique, simultnncous data from the Apolto IS and 16 surface magnetometers 
are used, requiring no orbital magnetic field data. This method has the advantages 
of using only the higher resolution surface magnetometers, and is insensitive to 
instrument otTsets. The lunar ionosphere in the geomagnetic tail is modeled, and 
results are compared with those of Apollo subsatcilite measurements (Kusscll vt 
ai, I974a,b). Values of lunar magnetic permeability and iron abundance are 
calculated taking into account effects of the diamagnetic ionosphere. 


Lunar Electrical Conductivity and Temperature 

Lun.'ir electrical conductivity is calculated using measurements of global eddy 
current fields induced by changes in the magnetic field external to the moon. The 
time dependence of the induced field response is a function of the electrical 
conductivity distribution in the lunar interior. Simultaneous measurements of the 
transient driving field and the lunar response field, by Explorer 35 and Apollo 
surface magnetometers, allow calculations of the conductivity. Detailed descrip- 
tions of the ApoHo and Explorer 35 instruments are reported in Dyal et al. (1970) 
and Senett et ai (1967). 


Analytical technique 

Data have been selected from measurements obtained in the lobes of the 
geomagnetic tail, during times when there is no indication of plasma effects 
(Anderson, 1965). The individual data sets are also subject to the following data 
selection criteria: (1) the magnitude of the field external to the moon (measured by 
the Explorer 35 magnetometer) is required to be at all times greater than 8 
gammas; (2) the external field is directed approximately along the :un-earth line; 
(3) main qualitative features of each event are required to appear in both surface 
and orbital data in all three vector coordinate axes to minimize use of data with 
large field gradients between the two magnetometers; and (4) no plasma data are 
measured above the solar-wind spectrometer instrument threshold. Criteria (1) 
and (2) a-e used to insure that no neutral sheet crossings are used in the analysis. 
Therefore, we have chosen data for times when the lunar response can be 
modeled by that of a conducting sphere in a vacuum. Details of the vacuum theory 
are reported in Dyal et al. (1972, 1974). 

Since the field equations governing the electromagnetic response of a sphere in 
a vacuum are linear, the sum of any two solutions is itself a solution. We take 
advantage of this principle by superimposing many driving functions (measured 
by the Explorer 35 orbiting magnetometer) to form an aggregate data event, and 
comparing this with the aggregate formed by superposition of the corresponding 
response functions (measured by an Apollo surface magnetometer). These 
aggregates have much higher signal-to-noise characteristics than do the individual 
events since errors such as those due to magnetometer digitization and to geotail 
gradients will be gaussian and tend to average out in the event-addition process. 
The response to the driving field aggregate is calculated for an assumed electrical 
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conductivity profile a(r), comptired to the measured time scries response (Apollo 
magnetometer data aggregate) and then reiterated by adjusting <7(r) until the error 
between the calculated response and measured response is minimized. 

Electrical condiiclivity results 

The eighteen data events -sed in our analysis have been linearly superimposed 
using a 20-inch IMLAC programmable display system in a real-time interactive 
mode with the Ames IBM 360/67 computer. For an assumed conductivity profile 
the IMLAC display is used to compare the calculated response with the measured 
Apollo aggregate. The conductivity profile is then iteratively adjusted until the 
calculated response function matches the measured Apollo response. Transient 
aggregates for the radial and two tangential axes are shown in Pigs. 1,2, and 3. The 
“computer response field" in each of these figures is calculated f*'om the 
conductivity profile displayed in Fig. 4. The figure insert shows a family tii profiles 
which also fit the data and give bounds on lunar conductivity. One of the profiles 
shown in the insert is the homogeneous-moon prefile, which gi'.es a value of 
7 X 10*^ for the average bulk conductivity of the moon. 

In our analysis we assume that the electrical conductivity increases monotoni- 
cally with depth from 10 * mho/m at the surface (Dyal and Parkin, 1971b) and that 


AGGREGATE TRANSIENT RESPONSE (RAOIAL-X) 



Fig. t. Aggregate transient response tradial-x). The aggregate is the sum of radial 
(ALSEP x-axis) magnetic field components for eighteen events in the geomagnetic tail 
lobes, selected to minimize plasma elTects. llie external field is a sum of events measured 
by lunar orbiting Explorer 35 mugneiomeier; measured response field is the surface 
magnetic field aggrcgaie of the ^ame events measured by the .Apollo 12 magnetometer. 
Computed response field is a iheorelicai response calculated for a sphere of electrical 
conductivity <r(r> shown in Fig. 4. 
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AGGREGATE TRANSIENT RESPONSE (TANGENTIAL- Y) 



Fig. 2. Aggregate transient response (langcniial-y). Aggregates are the external, com- 
puted, and measured response fields for the tangential (eastward) ALSEP V-axis 
corresponding to the radial A'-axis data illustrated in Fig. t. 


aggregate transient RESPONSE (TANGENTIAL -2) 



Fig. 3. Aggregate transient response (tangential-;). Aggregates are the external, com- 
puted, and measured response fields for the tangential (northward) ALSEP Z-axis 
corresponding to the radial A'-axis data illustrated in Fig. 1. 
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LUNAR ELECTRICAL CONDUCTIVITY 
DEPTH, km 



Fi{!. 4. Lunar electrical conductivity profile which gives “best-fit** response to the 
uggrcgale event shown in Figs. 1, 2, ond 3. Insert shows other conductivity profiles which 
give approximate bounds on the range of allowed solutions. 


it is a continuous function from the surface to the center of the moon. (The 
conductivity is shown for alt radii in Figs. 4 and S since the analysis requires that 
the conductivity be defined throughout the entire lunar sphere.) The uncertainty in 
the electrical conductivity is reflected in the spread of allowable profiles shown in 
Fig. 4 insett. These profiles are the simplest forms that we consider appropriate 
for our data set: the conductivity below 0.3R„, however, can vary by orders of 
magnitude without being discriminated by our technique. These uncertainties 
cannot be delineated by error b.irs associated with individual depths because 
errors in the analysis are reflected by the entire conductivity function. The 
uncertainty in the profile arises from several factors: (I) the nonuniqueness of 


LUN4H ELECTRICAL CONDUCTIVITY 
OCPTH, Iwn 

. 1500 1000 500 0 

iQ.tf— T , , 



R/Ry 


Fig. 3. Electrical conductivity profile illustrating the maximum size of a hypothetical 
highly conducting core of conductivity 7 x 10 ' mho/nu <s0.57K,. 
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profile* dctcrniinaticm as discussed by Backus and Gilbert (1970), rhillips (1972), 
and Hobbs (1973); (2) the penetration depth allowed by ti e lenelh of the data set 
used and its frcetiiency content; (3) the frequency response limitations of the 
Explorer 35 and Apollo magnetometers; (4) nonhornogcneitics in the external field 
over the dimensions of the moon; and (5) Instrumental errors in the measured 
fields. The conductivity profiles shown in Fig. 4 arc consistent with previous 
transient-analysis results (Dyal and J’arkin, 1971a; Dyal et al„ 1974). 

An important result of this analysis is that the lunar conductivity rises rapidly 
to 10*' mho/m within 250-kin depth and remains relatively constant to 800-kin 
depth. Wc have also investigated the possibility that a partially molten core exists 
in the moon at a depth of approximately 800 km, as indicated by seismic results 
(Nakamura c( ai, 1974). We have examined the limits that electrical conductivity 
analysis places upon a partially molten core, using an order-of-magnitude increase 
in conductivity to represent a phase change (Khitarov et til., 1970; Fresnall ct ai, 
1972), then determining where such a highly conducting core could exist consis- 
tent with magnetic field data. The upper limit on the size of such a core is 
determined to b^; 0.57 lunar radius, close to the value 0.55fin,<»n determined from 
!c measurements (Nakantura et at,, 1974). Analysis of longer transient-event 
vgG' -etdes will further define the size limit on a highly conducting core in the 
moon. 

Inferred lunar temperature 

For minerals which are probable constituents of the lunar interior, the 
electrical conductivity can be expressed as a function of temperature, and 
therefore a thermal profile of the lunar interior can be inferred from the 
conductivity profile in Fig. 3 for an assumed material composition. Many 
investigators have published laboratory results relating conductivity to tempera- 
ture for materials which are geochemical candidates for the lunar interior le.g.. 
England et at., 1968; Duba et at., 1972; Schwerer et a!., 1972; Duba and Ringwood, 
1973; Olhocft et at., 1973). 

Recently Duba et at. (1974) have measured conductivity of olivine as a 
function of temperature under controlled oxygen fugacity and have found the 
measurements to be essentially pressure independent up to 8 kbur. These meas- 
urements for olivine have been used to convert our conductivity proiile of Fig. 4 
to a temperature profile (Fig. 6). The Fig. 6 insert shows a family of temperature 
profiles which reflect the error limits of our conductivity calculations. We see 
from Fig. 6 that the iron Curie isotherm (--780"C) should be reached within a 
depth of 200 km, allowing permanently magnetized material to exist only at 
shallow depths in the moon. The results imply a high temperature gradient to 
250-km depth (averaging 6‘'K/km), and rather uniform temperature from 250 to 
800-km depth. Similar conclusions have been reached by Kuckes (1974) using 
harmonic analysis of magnetometer data with the interpretation that convection is 
an important process in the present day moon. Also Hanks and Anderson (1972) 
have theoretically calculated a thin thermal boundary layer and high temperature 
interior for the moon, consistent with our results. 
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LUNAR TEMPERATURE 



F>e< 6- Temperature profiles of the lunar interior, determined from the conductivity 
profiles in Fig. 4. A lunar interior composition of olivine is assumed: laboratory data for 
olivine from Duba el of. (1974) are used in the calculation. 


Lunar Magnetic Permeability and Iron Abundance 

Relative magnetic permeability of tlie moon can be calculared using data 
measured during times when the global magnetization field is the dominant 
induced lunar field. This type of induction dominates when the moon is in 
magnetically quiet regions of the geomagnetic tail, where the ambient external 
field is essentially constant and therefore eddy current induction is minimal. 


Analytical ft'c/iniqui? 

In previous reports (Parkin et ai, 1973, 1974) we calculated permeability p. 
using simultaneous measurements of the external field H by the lunar orbiting 
Explorer 35 magnetometer and the total surface field B using the Apollo 12 or 15 
lunar surface magnetometer (LSM). In the present analytical n.ethod, simultane- 
ous magnetic field data from Apollo 15 and 16 LSM's are used, requiring no orbital 
magnetometer data. This method has the advantages of ( 1) higher resolution of the 
LSM*s and (2) results are not sensitive to instrument offsets, since calculations 
involve field differences, and offsets are canceled out of the analytical equations. 
Also, data are selected from deep in geomagnetic tail lobes where plasma eflfects 
are minimized, and lO-min averages of steady data are used so that eddy current 
induction effects are negligible. 

Our present method involves using sets of field measurements D made at two 
times, each simultaneously at mo Apollo sites (by the Apollo 15 and 16 LSM’s). A 
description of the theoretical development is given in Appendix A. In this method 
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a system of twelve equations is used plus a transformation matrix relating Helds 
measured at the two ::uifnce sites, yielding a total of 13 equations with 10 
unknowns. These equations are solved for the lunar bulk relative magnetic 
permeability p.i,, expressed as follows: 


AJf “ OtiABi, 

0) 



(2) 

AB}, — 0::A/iiy — ai)AB u 

fljiAB It 

(3) 



where subscripts I and 2 denote Apollo sites, e.g. site I can be designated Apollo 
15 and site 2, Apollo 16; subscripts .v, y, z, denote vector components in the 
ALSEP coordinate system with oiigin at either site, where x is directed radially 
outward from the lunar surface, and y and f arc tangent to the surface, directed 
eastward and northward, respectively: «</ is an element in the transformation 
matrix from site I to site 2. Each Held-ditTcrence term AB,, => BuiU) 

denotes a dilTcrcncc in a field component measured at the same site at two 
different times when field values are different. We use the generalized notation 
fii ~ Ai/Aj, where Ai and A; are defined as numerator and denominator, respec- 
tively, of cither Eq. (I), (2), or (3). 

S 

Rela . s masnetic permeability results 

To solve for bulk relative permeability of the moon we perform a regression 
analysis on the generalized field-dilTerence parameters Ai and As, shown in Fig. 7. 
This curve has been constructed using over 2000 simultaneous Apollo 15 and 
Apollo 16 magnetometer data sets, measured during five orbits of the moon 
through the geomagnetic tail. These data have been carefully selected to eliminate 
data measured in the plasma sheet or contaminated by other induction modes such 
as eddy current induction (see Parkin et ai, 1974). 

Effects of plasma diamagnetism and confinement are minimized by eliminating 
data points for which the magnitude of the external magnetizing field H < 
7 X 10“’ Oe. Since poloidal eddy current induction is dependent upon time rate of 
change of II, inclusion of poloidal fields is minimized by averaging over 10-min 
intervals during which Apollo 15 and 16 data peak-to-peak variations are small. 
Then pairs of data points are selected from diff erent times, using criteria to insure 
that the denominators in Eqs. ( 1) to (3) are non-zero. Finally, a regression analysis 
is performed on the set of ordered pairs (Ai, A;) plotted in Fig. 7. The least-squares 
best estimate of the slope, calculated using the method of York (1966), is 

= 1.008 ±0.005. 

Lunar magnetic permeability results adjusted for ionospheric effects 

We have calculated a value for relative magnetic permeability of the moon. To 
determine the absolute lun:tr permeability and iron abundance in the moon, we 
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LUNAR PERMEABILITY ANALYSIS 



Fig. 7. Lunnr permeability analysis, using all three vector axes of simultaneous Apollo 
IS and 16 surface magnetometer dam for live lunations in the geomagnetic tail lobes. The 
least-squares slope gives a value for the bulk lunar relative permeability: 1.008 s: 0.005, 
The parameters and iS: are defined in the text. 


consider the magnetic permeability of the environment exterior to the moon, in 
particular, the lunar ionosphere in the geomagnetic tail. 

The lunar atmosphere (see Johnson, 1971; Hodges et ai, 1974) is the source of 
the ionosphere. The lunar atmosphere is an exosphere whose sources are 
neutralized solar-wind ions and neutral atoms outgassed from the moon. The 
species of solar-wind origin are ’‘Ar. '“Ne, ‘‘He. H:, and H. Hydrogen and helium 
escape thermally with lifetimes ranging from 10' to 10’ sec while ’"Ar and ‘“Ne 
have much longer lifetimes which are controlled by photoionization. 

Previous work on the lunar ionosphere has been confined to time periods when 
the moon is in the solar wind (Manka, 1972; Vondrak and Freeman, 1974). During 
the four-day period when the moon is in the geomagnetic tail, the solar wind is no 
longer a primary source of neutral atmosphere. He, H:, and H thermally escape 
within several hours and do not contribute significantly to the lunar ionosphere. 
*®Ne, “Ar, and ’"Ar are the main atmospheric <- ^nstituents due to their long 
residence times. The ionization of the atmosphere .. he solar wind is dominated 
by charge excliange (Weil and Barasch, 1963), whereas in the geomagnetic tail the 
principal process is photoiomzaiion by ultraviolet radiation from the sun. In 
addition, acceleration by the motional electric held has been considered to be the 
principal loss mechanism for ionospheric particles when the moon is in the solar 
wind. In the geomagnetic tail, plasma conditions are different; thermal escape, 
electric field acceleration, and ballistic collision with the surface are possible loss 
mechanisms. 

The photoionization of the sunlit lunar atmosphere is determined by the solar 
spectrum and the photoproduction cross sections of each species. In this 
photoionization process the photon energy is used to ionize the neutral atoms and 
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the excess energy is converted to pliotoelectron kinetic energy. By considering 
the peak photon flux and cross sections, we estimate the photodectron energy to 
be about 20 cV for both Ne and Ar. We assume that the ions remain at the 
effective temperature of the neutra) atmosphere which is in thermal equilibrium 
with the surface at 300“K. The thermal speeds of the high energy electrons arc in 
excess of the escape velocity (2.4 x lO’ cm/sec). However, the electrons arc also 
electrically attracted to the ions which arc gravitationally bound to the moon at 
the low temperature of 3()0°K, To consider these competing cifccts we consider 
the equations governing the behavior of both the ions and electrons of each 


ionized species: 


1 d«*.. -g(r)m' eE 
It* dz " tcT* kT* 


(4) 


1 dn,^ -g(r)m, eE 
n.dz kT, ~kT, 

where n*, in*, T*, and «„ in„ T, are the density, mass, and iemperature of the 
ions and electrons, respectively; z is altitude above the lunar surface; g(r) is the 
lunar gravitational acceleration; e is charge on the electron; and k is Boltzmann’s 
constant. The electric field E is the Rosscland field, which is responsible for 
maintaining local charge neutrality. Since g(r) takes into account the finite mass 
of the lunar atmosphere, solution of Eqs. (4) and (5) for the density is valid far 
from the moon as well as near the lunar surface. A scale height can be calculated 
from this distribution which depends, for the model described above, :>n the ion 
mass and the electron temperature and is about 1 100 km for both neon and argon. 
Using an e.\prcssion given by Johnson (1971) for the mean residence time t, we 
can now estimate the loss rate of the ionosphere by thermal escape. We also 
calculate the production rate for the ionospheric constituents. Using this informa- 
tion wc then calculate the time derivative of the ion density: 

~«p-n/T (6) 

where p is the production rate and n It is the loss rate. At equiliarium n = pr. 

Results indicate a characteristic ion density of approximately 10 cm”’ between 
the surface and 100-km altitude. The energy density of this ionospheric plasma is 
dominated by the energetic electrons, and the geomagnetic tail field is about 10 
gammas. The ionospheric plasma diamagnetic permeability p, is derived from 
B = pH = H +47 tM, where M = - nfcTB/B’: 

p, =(l + .a/2)-' (7) 

where /3 -EnnkTlB^. Using ion density n - 10 ions/cm’, plasma temperature 
Ts T, = 1.5x 10’®K and B = 10 gammas, we calculate the ionospheric relative 
magnetic permeability to be p, = 0.8. This theoretical diamagnetic permeability 
for the ionosphere compares favorably with the experimental value obtained by 
simultaneously considering our lunar relative permeability result and the experi- 
mental lunar induced dipole moment determined by Russell et al. (1974a). A 
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two-layer (ionosphere-moon) muilcl is nsec! to calculate permeability of tlte moon 
am! the ionospliere relative to free space (see Parkin ct al„ 1974). The ionosphere 
diamap.netic permeability is calctihited to be fh =0.761o li and the paramagnetic 
permeability of the moon, adjusted for ionospheric citects, is fi = 1.012 This 
rcstill is in general agreement with previous measurements using Explorer 35 and 
Apollo 12 nmgnetometcr data reported by Parkin et at. (1974). 


Iron abimdaiice 

Using the value of global lunar magnetic permeability, we can determine free 
iron and total iron abundances. The free and total iron values are also constrained 
by lunar density and moment of inertia, and are functions of thermal and 
compositional models of the lunar interior. 

The lunar bulk permeability n « I.OI2-SI«vi is too high to be accounted for by 
any paramagnetic mineral which is a likely constituent of the lunar interior, 
implying that some material inside the moon must be in the ferromagnetic slate. 
Assuming the ferromagnetic material is free iron of nonintcracting multidomain 
grains, the lunar free iron abundance can be determined using a thermal model of 
the lunar interior. The thermal proftlc is approximated by a two-layer model with 
the boundary located at the iron Curie point isotherm. Figures 8 and 9 show free 
iron abundance {q) and total iron abundance (Q) related to (/x — I) and thermal 
profile. Q is shown for two compositional models; we assume the moon is 
composed of a homogeneous mineral (olivine or orthopyroxene) of uniform 
density 3.34 g/cm^ with free iron grains disbursed uniformly throughout the 
sphere. All assumptions and the theoretical formulation leading to Figs, 8 and 9 
are outlined in Appendix B. Using the temperature profiles shown in Fig. 6, we 
find that the iron Curie isotherm radius Rr should be in the range A 0.88, where 
A = RjRm. In our calculations we use A =0.88. 

Without adjusting our measured bulk permeability for ionospheric effects we 
determine free iron abundance to be 1.6± 1.0wt.9£. This corresponds to a total 
iron abundance of 6.7£0,6wt.% for the free iron/olivine model and 13.2 c: 
0,4 wt.% for the free iron/orthopyroxene model. Using fi = l.0l2!,n«,i, which does 
account for the lunar ionosphere, we determine free iron abundance to be 

2.5 m:? wt.%. The free iron abundance corresponds to total iron abundance of 
6.0 ± 1.0 wt.% for the free iron/olivine lunar model, or 12.5 ± 1.0 wt.% for the free 
iron/or*hopyroxenc model. If we assume the lunar composition to be one or a 
combination of these minerals, the total iron abundance will be between S.O and 

13.5 wt.%. 


Summary and Conclusions 
Lunar electrical conductivity and temperature 

(1) Conductivity results presented in this paper have been determined using a 
new transient-superposition technique which increases the data signal-to-noise 
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Fig. 8. Free and lota) iron abundances in the moon as a function of the built magnetic 
permeability n and temperature profile of the moon, for a free iron/olivine compositional 

model. 

ratio and uses both radial and tangential response data in the geomagnetic tail 
lobes. 

(2) The average bulk electrical conductivity of the moon is calculated to be 
7 X lO"* mho/m. 

(3) The calculated radially varying conductivity profile fsee Fig. 4) rises 
rapidly with depth to ~ 10 ' mho/m within the first 250 km, then remains relatively 
constant to 800-km depth. 

(4) The limiting maximum size of a highly conducting core (3=7 x 10"’ mho/m> 
is calculated to be ^O.SIR/Rm^n: that is, our present resolution allows us to probe 
to depths where moonquakc foci have been located by seismic studies. 

(5) A conversion of electrical conductivity to temperature using the data of 
Duba et a!. (1974) for olivine yields a thermal profile which is relatively hot, 
implying that the Curie point is within 2G0 km of the lunar surface. 

Lunar magnetic pemie<ibii//y and iron abundance 

(1) Magnetic permeability results have been obtained using simultaneous data 
from the Apollo 15 and 16 surface magnetometers and requiring no orbital 
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IRON ABUNDANCE, FREE I RON /ORTHOPYROXENE MOON 
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Fig. 9. Free and total iron abundances in the moon for a free iron/orthopyroxene lunar 

model. 


magnetic field data. Both radial and tangential magnetic data are used in the 
analysis. 

(2) The bulk magnetic permeability of the moon relative to its environment is 
calculated directly from magnetometer data to be fix = 1.008^ 0,005. 

(3) The lunar ionosphere in the geomagnetic tail is modeled, and theoretical 
calculations yield an ionospheric permeability of 0.8, a value consistent with 
surface and orbital magnetometer measurements. 

(4) The magnetic permeability of the moon, adjusted for effects of its 
diamagnetic environment, is calculated to be ^ = 1. 01 2 -O IH1N. 

(5) Free iron abundance is determined from our permeability and temperature 
calculations to be 2.5r; ?wt.%. 

(6) Total iron abundance is found for two compositional models of the moon: 
free iron/olivine model. 6.0±1.0wt.^; free iron/orthopyro,\ene model, 12.5 ± 
l.0wt,%. Assuming the moon is composed of one or a combination of these 
minerals, the overall iron abundance will be between 5.0 and 13.5 wt.%. 
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Appendix A 

It has been shown by Parkin el of. (1973) that the field on the surface of a spherically symmetric 
two-layer magnetically permeable sphere can be expressed as; 

B = H.(l 2F)i + JLd - F)p + H.(l - F)f 

where 

p (2r) + n(M.-l)-A*(T)-D(2M.-M) 

“{2ij + l)((a. + 2)-2A’(n-l)()i.-l) 

Here H is the field e.Ylernal to the sphere at large distances; t) - and are relative magnetic 

permeability of the shell and core respectively (permeability of free space n„ = 1); A = R;/R, where Ki 
is the radius of the boundary between the two permeable regions and R , is the radius of the sphere. .As 
applied to the lunar sphere Eq. d) is e.xpresscd in the ALSEP coordinate system which has its origin 
on the lunar surface. The .t axis is directed radially outward from the surface; the y and ; axes are 
tangential to the surface, directed eastward and northward, respectively. 

In order to measure the dipole induced in Ihe moon by the external field H w-c subtract two values 
of B in Eq. d) for which 11 is different: AU = B-B'. Forming this difference at two dillerent points on 
the moon for the two values of the e.xternal field we have in component form; 

Afl„ = d + 2F)AH.. 

AB., = (1-F)AH„ (3) 


(I) 

( 2 : 
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where f 1. 2 for the two ditfcrenl mai;nctomctcr sites. The two Helds Hi can be rehited by the 
Irunsfotncition 

H.-lAin. 

II.-lAl 'H, 

where |Al is iMc tinnsformnlion matrix relatinp the ALSHP coordinate systems at the two sites. 
Equations (.t) and (4) can be solved sinmltancously eliminating II, to give 


ABi, “ 


/I A B 1 1 , AB, 1 , AB I, I 

(1 + 2F)[^o„ ^ ,pj+ 0,5 ^ j I _ 

(1 - F)|^o„ -j + 2D (I - F) " (I - F jj 

n ABu , , AB,. , ^ AB„ 1 

(1 - Dj^o,, ^ (1 _ (1 _ f , J 


(5) 


where a,i nrc the elements of |A|. From Eqs. (5) we have three expressions involving the induced 
magnetic moment l/Ki*F which are: 


with 


■ I S3 I 

* OtvAi^ii *1' 

OjiABti 

/i* 




A0}p — ~ QjySBit 

ABir “ fljiAHiy •” OjiABij 


t + ?F 

•*““T:f 


(6) 


(7) 


since Eqs. (5) could have been expressed with AD, as the dependent variable, Eqs. (6) can also be 
written interchanging I and 2 subscripts on the AB quantities and using the matrix elements of |A| ' 
instead of those from |A|. 


Appendix B 

Limits imposed on (i. of Eqs. (6) in Appendix A by the LSM data can be used to calculate the lunar 
iron abundance for suitable lunar compositional and thermal models. Two of the models are described, 
including magnetic and other geophysical constraints. In both cases the lunar interior is modeled by a 
sphere of homogeneous composition. Free iron of multidomain noninleractiri; grains is assumed to be 
uniformly distributed throughout a paramagnetic mineral. The paramagnetic component in one case is 
olivine Jy Fe.SiO, • (1 - y)Mg:Si 04 | and in the other case is orthopyroxene (y FeSiO, (I - y)MgSiO,l. 
The free iron is ferromagnetic in the regions where the temperature T is less than the iron Curie 
temperature T, and it is paramagnetic w'hcre T > T,. Therefore, each model is a two layer permeable 
sphere where #t, and /i, arc the relative magnetic permeability of the spherically symmetric shei! and 
core respectively; B. is the core-shell boundary. I Equation (2) of Appendix A rehaics these variables to 
the measured magnetic moment HKi'F.I We consider the magnetic contributions from both free and 
combined iron in both the shell and core of the model with p., : ■= I -t- 4rrK,., where 

K,-K,(y.T,)-hKV(q) 

K.-Kp(y.T,)-i-K^(q) 

and Kf is the paramagnetic susceptibility of the olivine or orthopyroxene, y is the mole fraction of 
ferrosilicate in the mineral, T, is the uniform temperature of the shell, T- is the uniform temperature of 
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the core, K, is the apparent ferrumaynctic susceptibility of free iron in cinu/cm', K f Is me apparent 
paromafiieiic snseeptibility of free iron (above the Curie tcmperatcre). tmd q is the mass fraction of 
free iron in the moon. 

The measured susceptibility of the free iron K, is an apparent value which differs from the intrinsic 
susccptibiliiy of iron K because of self-demapnctbaiiun of the iron grains and the fraction of iron in 
the moon. 1'lie apparent ami intrinsic susceptibility of the free iron are related (see Niiga(a, )96l) by 




qp K 

Pr\ + NK 


( 2 ) 


where p is (he lunar dcnsi(y (assumed uniform throughoul the moon), and pi is (he densi(y of iron. An 
analogous expression relates fs|. (he apparent paramagnetic susceptibility, and K', lire intrinsic 
paramagnetic susceptibility of free iron (T > T. ), 

Nagata e( al. (I957) found the susceptibility of olivine to be 2.4y x lO 'emu/motc at room 
leinpcraluro. Using this c.xprcssion, the Curie law temperature dependence, and the following 
cnipiricnl equation for olivine density (Dana, 1966) 


pp(y) •= 0.94y + 3.26 g/cm’, (3) 

we obtain for olivine 

K, O' 0. 14 ^ gffrilf W) 

Similurty Akimoto ft al. (1958) give the susceptibility for pyroxenes; l.ly x lO'^emu/mote. From Deer 
tl al (1962) we obtain the following empirical expression for pyroxene density; 


pp(y) = 0.71> +3.15 g/cm'. 


(53 


Again, combining the expressions from Akimoto et al. and Deer et al. with the Curie law temperature 
dependence, we obtain for pyroxene 


( 6 ) 


The lunar moment of inertia is approximately that of a sphere of uniform density (Ringwood and 
Essene, 1970). Choosing a uniform density for our lunar model we can write 


p pr Pf(y) 


(7) 


The free iron abundance q can be determined as a function of p for the orthopyro.xcne/free iron model 
by simultaneously solving F.qs. (2) and (7) of Appendix A and Eq... (I),(2). (5), (61, and (7) above. The 
olivine/free iron model solution is determined by solving Eqs, (2) and (7) of Appendix A and Eqs. (I), 
(2), <3), (4), and (7) above. Results for iron abundances as a function of lunar magnetic permeability are 
given in Figs. 8 and 9. <p, the fractional lunar mass due to both chemically uncombined and combined 
iron, is constrained by q and the lunar density. For the orthopyroxene/free iron model 


® “'*‘''32y + 100^* 

(8) 

and for the olivine/free iron model 


- . Il7v /. 
<?“'>-"64y + 140^'-‘'>- 

(9) 


In Eqs. (8) and (9) y can be eliminated by substitution of F,qs. (3). (5), and (7); therefore, the total iron 
and free iron abundances are directly related. Figures 8 and 9 show the results of these culcuiutions 
relating q, Q and p - I. 
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ConMunis u»cd in llic calculntions for both models ore; 


uniform dciisliy of moon 

density of free iron 

denuittneliMtiuii f;icior of iron prains 

initial inti insie ferromagnetic suscci'libiliiy of iron 

intrinsic paramapnclic susceptibility of iron 


thermid models 


3.34 g/cm’ 

7.8Sg/cni’ 

N=3.5 (Nogatu, 1961) 

K» I2cnui/cm‘ (Bozorth, 1951) 
K"«’2,2x |0 *cmu/cm’ 

(Tebble nnd Craik. 1969) 
T,“900‘K, T,= I700'K. 

A » 0.98. 0.96, 0.94, 0.92 
T, ■= 800’K. T, « I600’'K. 

A » 0.88, 0.84 
T, « 700’K, Tt “ 1400"K. 

A-' 0.80 

r,« 600 "K, T,*= IOOO‘K, 

A -0.70 


s 
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STRUCTURE OP THE LTOIAR nrTERIOR FROM MACI®TIC FIELD MSASUREI-ffiriTS, 

P. Dyal*, C.W. Parkin**, and W.D. Dally***. 

*NASA-Ames Research Center, Moffett Field, CA 9^035 
**Univ. of Santa Clara, Sarita Clara, CA 95053 s 

***Eyring Research Institute, Provo, UT 8J46OI 

The electrical conductivity and magnetic permeability of the lunar 
Interior have been determined from measurements by a total of six lunar sur- 
face and orbiting magnetometers. From these results*, characteristics of 
lunar internal structure are inferred. 

A new technique has been applied to conductivity analysis in which 
simultaneous data is used from a network of three instru-ments: The Apollo 

15 lunar surface magnetometer (ISM , the Apollo I6 LSM, and the Apollo I6 
subsatellite magnetometer, which provides coverage around the entire global 
circumference. Also, a two-instrument technique has been applied which uses 
Apollo 12 LSM data and simultaneous data from the Goddard magnetometer aboard 
Explorer 35* All measurements are made when the moon is located in 
high-latitude regions of the geomagnetic tail where plasma effects in the 
lun a^* environment are minimal. Individual magnetic events are superimposed 
to obtain a single large transient for analysis. Also, dui'ing examination of 
five years of data one exceptionally large, well-behaved transient was 
recorded when the moon was in the geomagnetic tail. This single event 
allowed substantial improvement in resolution and sounding depth for conduc- 
tivity analysis. 

Although the t;»o types of conductivity analj'^sis are analytically 
different and use measurements obtained over a five-year period from six 
different magnetometers, results are in surprisingly close agreement. The 
conductivity profile result is shown in the accompanying flgui-e. A striking 
feature is the abrupt transition near 300 km depth where a knee occurs in the 
conductivity profile. The conductivity increases rapidly from the surface to 
about 4xlO“J ^os/m at 300 km depth. At greater depths the conductivity 
Increases more slowly to about 2x10”^ rahos/m at 8C0 km depth. 

This conductivity transition at 300 km depth, the location of which 
corresponds closely to that of the seismic velocity change reported by 
Nakamura et al. (1974), strongly implies a structural or compositional 
change at that depth. 


In addition, the lunar magnetic permeability determined from magnetome- 
ter measurements has been used to place limits oh a possible higlily ccnduut- 
Ing core in the moon. In the permeability analysis the moon is represented 
,by a three-layer magnetic model; an outer shell of temperature (T) belo'w the 
Curie point (T ), whose permeability u is dominated by ferromagnetic free 
iron; an intermediate shell of T > T where permeability » that of free 
space; end a V.ir'nly conduaoing cere ?c » 10"" mhor/m.' mcdslad by u = 0. 

Tills core effectively excluceo c-xesrna-l maj.ieeic fields ever time lengths of 
days and therefore acts as a strongly diamagnetic region (u -»0). 



A theoretical analysis naa oeen carrxea out to relate tne inaucea mats- 
netic dipole moment to the core size. The induced dipo' moment has been 
detemlned from simultaneous Apollo 12 and Ames Explorer measurements to be 
2.1x10^^ ± l.OxlO-'*'^ gauss cm^ (Parkin et 1974) end from simultaneous 
Apollo 15 and 16 measurements to be 1.4x10^^ ± 0.9x10^® gauss cm3 (Dyal et al.j 
1975). 

Results show that the core size is a function of the depth of the Curie 
isotherm and lunar composition. To calculate an upper limit on core size, we 
assume that the compositional mineral of the moon has no iron silicate and 
take an extreme depth of 350 km for the Curie point (Toksoz et al. , 1974) so 
that free iron content is maximized. Then the maximum allcv.’’able core size is 
470 km for an olivine moon and 600 ka for an orthopjrroxene moon. However, 
both magnetic permeability and electrical conductivity analyses are consistent 
with the absence of a highly conducting lunar core. 
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LUJIAK r.Lb’CTRlCAL C0:;DU9TIV15’V PROFILE FROM 
MAGEKTOMB'JER LPJIV.’CRK 

W. P. P''.l.ly (Eyrinc Kenearch Institute, Provo, 

UT 8!;60l) 

P. Pyal (K.^.SA-.Awc Research Center, Moffett Field, 

CA 9^1035) 

C. V/. Parkin (University of Santa Clara, Santa 

Clara, CA 95053) 

I'hc electrical conductivity of tlic lunar 
interior has been dcterniinca v;ith o. new analyti- 
cal technique in \,'hich sii;:Ultai^eous data ca'c used 
from a network of the lunar iriaGneter.eters; the 
Apollo 15 loiiar surface maf-^ietomster (LSK), the 
Apollo l6 IjSM, and the Apollo l6 subcatellite 
maEnetometcr. From tliis network mcasui'cnents of 
induced poloidal fields con be made over the 
lunar sphere, which are therefore rcpre.sentative 
of the v.’h ole* moon electrical cui'rent distribU' >.on. 
Magnetic field r.ieasurc.n-.ents v,'ere made dui'ing times 
when the moon was at hirh latitudes in the geo- 
magiietic tail lobes. Since in the lobes the 
plasma effects are mdnimal, the moon has been 
modeled as a sphere in a vacuum. In the analysis 
many individual magnetic events are superimposed 
to obtain a single large traJisient with much 
better signal- to-ncise characteristics than single 
events. The conductivity profile obtained from 
this 3-mag.neto.meter technique is in close agree- 
ment \rith previous results from our 2-magnctom.eter 
analysis using the Apollo 12 LSM and Explorer 35 
magnetometers. The conductivity increases rapidly 
from the surfe''.i to about Uxl0“3 nhos/m at 3C0 km 
depth. At greater depths tiie conductivity 
increases more slowly to about 2x10“'^ mhos/m at 
800 km depth. A striking feature of the profile 
is the abrupt transition near 3C0 Lm depth vfnere 
a knee occurs in the conductivity profile. This 
conductivity transition, the location of v/hioh 
corresponds closely to that of the lunar seisniic 
velocity change reported by IJakajr.ura ct al. (I97^i 
strongly implies structural or compositional 
changes at that depth. 


